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Abstract 

We investigate the nonhnear rheological behavior of colloidal suspensions of Lapo- 
nite, a synthetic clay, driven by a steady and homogeneous shear strain. We show 
that the external drive leads to a drastic slowing down of the aging dynamics or even, 
in some cases, in the rejuvenation of the system. Under shear, after a surprisingly 
long time, the spontaneous aging process observed at rest is suppressed. The system 
then reaches a non-equilibrium stationary state, characterised by a complex viscosity 
depending on the applied shear rate. In addition, the glass exhibits a non-Newtonian 
shear-thinning behavior. These rheological behaviors confirm recent numerical and 
theoretical predictions. 

1 Introduction 

Complex fluids, such as foams, emulsions, pastes and slurries, are knovv^n to exhibit a rich 
phenomenology of nonlinear rheological behaviors (Larson, 1999). In particular, similar low 
frequency shear rheology is observed for a v^ide range of these soft materials (see Fielding et 
al., 2000 for a review). For many of these materials, the stress response a to a shear strain 
of constant rate 7 is described hj a = A + ^7" that leads to a non-newtonian viscosity 
(Barnes et al, 1989). This common rheological behavior was recently attributed to two 
properties of these materials, structural disorder and metastability, resulting in a "glassy" 
dynamics, where thermal motion alone is not sufficient to achieve complete structural 
relaxation (Sollich et al, 1997; Sollich, 1998). 

Due to long-range interactions, it turns out to be a complex task to predict the nonlin- 
ear rheological behavior, including a non-Newtonian viscosity, on the basis of the structure 
and interactions in these systems. Numerical simulations have allowed to predict the re- 
lationship between the material's microstructure and its macroscopic properties (Brady et 
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ai, 2000; Melrose et al, 1995; Melrose et al, 1999). Recently, a different approach to the 
rheology of these "soft glassy materials" has been developed, starting from a model of a 
glassy system instead of taking all the interactions - including colloidal and hydrodynamic 
interactions - into account. Because of both their practical and fundamental interest, the 
nonlinear rheology of such soft glassy systems subjected to an external drive has recently 
attracted particular attention. This was done either by extending approaches for the de- 
scription of glassy dynamics (Sollich et al., 1997; Cugliandolo et al, 1997; Sollich, 1998; 
Berthier et al., 2000; Barrat and Berthier, 2000; Fielding et al, 2000) or based on phe- 
nomenological approaches (Hebraud et al, 1998; Derec et al, 2001, Coussot et al, 2002). 
The most interesting situation occurs when the relaxation time scale of the system becomes 
of the same order of magnitude as the time scale introduced by the shearing process. This 
usually happens in glassy systems because the typical relaxation time grows together with 
the waiting time tyj, that is the time elapsed after a quench from the "liquid" into the 
"glassy" state. 

At rest, the properties of glassy systems thus depend on the waiting time: the system 
is said to age. Recent theoretical and numerical works, initially developped for supercooled 
liquids, have allowed for a first detailed description of the aging process (Bouchaud et al, 
1996; Kob and Barrat, 1997). Glassy systems exhibit two modes of relaxation : a fast mode 
corresponding to the fast motion of particles in the 'cages' constituted by their neighbors 
and a slow mode of relaxation resulting from the structural rearrangement of these 'cages'. 
The strongest experimental evidence for the applicability of such predictions (including 
classical mode-coupling approach) comes from colloidal glasses (Nelson and Allen, 1995; 
Pusey and Van Megen, 1987; Bonn et al, 1999; Abou et al, 2001). Because of the 
experimental advantages they present compared to structural glasses, where there are two 
coupled control parameters, temperature and density, these colloidal glasses have been 
studied extensively. 

In this paper, we investigate experimentally the response of a colloidal glass of Laponite 
to a steady shear strain of constant rate. The experimental materials and methods are 
described in Section II. Section III presents the rheological experiments. We discuss our 
experimental results and conclude in Section IV. 

2 Materials and methods 

The experiments were performed on Laponite RD, a synthetic clay manufactured by 
Laporte Industry. The chemical composition of the clay is as follows : Si02, 65.82%; 



2 



MgO, 30.15%; Na20, 3.20%; Li02, 0.83% and corresponds to the chemical formula Sis.oo 
(Mg5.45Lio.4o)H4024Nao.75. The particles of Laponite are colloidal disks of 25 nm diameter 
and 1 nm thickness, with a negative surface charge on both faces (Kroon et al, 1998). 
The clay powder was mixed in ultra-pure water and NaCl at different clay concentrations 
varying between 2% wt to 5% wt. The ionic strength of the suspension, denoted /, was 
subsequently adjusted by dissolution of NaCl in water, in the range 10~^ < I < 8.10^^ 
M. The pH value of the suspensions was fixed to pH = 10 by addition of NaOH, provid- 
ing chemically stable particles (Thompson and Butterworth, 1992; Mourchid and Levitz, 

1998) . The suspension was stirred vigorously during 15 minutes and then filtered through 
a Millipore Millex - AA 0.8 /im filter unit. This preparation procedure allows us to obtain 
a reproducible initial liquid state. The aging time = of the suspension is defined as 
the moment it passes through the filter. Within a time varying from a few minutes to 
a few hours for the different concentrations and ionic strengths considered here, a three 
order of magnitude increase in the suspension viscosity was observed. Being in a liquid 
state right after preparation, the suspension becomes more and more viscous and does not 
flow anymore when tumbling the recipient. Since the physical properties of the suspension 
depend on the time elapsed since preparation t^, the sample is said to age. The aging 
dynamics can be characterised by oscillatory shear experiments in order not to disturb 
the system. These measurements are performed using small oscillations at a frequency of 
1 Hz at an imposed stress of 0.1 Pa. They yield G' and G", the elastic and loss moduli 
respectively. The complex viscosity modulus rj* can be calculated from these quantities 
from 7]* = (G"^ + G"'^y/'^ /uj. It can be observed that the viscosity changes very rapidly, 
by 2 orders of magnitude over about 1 hour, as shown in Figure [H Thus, suspensions 
of Laponite become strongly visco-elastic in time, even at low volume fraction $ in par- 
ticles (of the order of $ ~ 0.01). The formation of a gel, evidenced by the existence of 
a fractal network, has been invoked in explaining the visco-elasticity (Pignon a/., 1997; 
Mourchid and Levitz, 1998). Recently, the structure and viscosity of Laponite suspensions 
at ionic strength / = 10~^ M were studied by using static light scattering (Bonn et al, 

1999) . Contrary to previous observations, no evidence for a fractal-like organisation of the 
particles could be found, provided the suspensions were flltered. Therefore, it has been 
proposed that Laponite suspensions form so-called repulsive colloidal glasses, stabilised by 
electrostatic repulsions. In such systems, the ergodicity is lost due to blocking of particle 
movements by the dense surrounding cages formed by their nearest neighbours. The col- 
loidal glass is obtained for very low volume fraction $ ~ 0.01 compared to those for usual 
spherical colloids, for which glasses are obtained above $ ~ 0.5 (Pusey and Van Megen, 
1987). To account for this difference, one has to consider the effective volume fraction of a 
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Figure 1: Complex viscosity as a function of the aging time ty, for a Laponite suspension 
at 1.5% wt, / = 7.10-3 M. 
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Figure 2: Effect of the addition of NaCl on the driven aging dynamics : viscosity as a 
function of the aging time tw for Laponite suspensions at 1.5% wt, / = 5.10"^ M (cross) 
and / = 7.10"^ M (circles), by applying a constant shear rate 7 = 500 s"^. 
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particle, which is done by adding to the particle radius the Debye length. To give an order 
of magnitude, it is estimated to 30 nm for a ionic strength / = 10""^ M. Recent experi- 
ments have shown that the location of the glass transition line in the (volume fraction / 
electrolyte concentration) phase diagram is consistent with this assumption (Levitz et al, 

2000) . 

Addition of salt to a suspension of Laponite disks reduces electrostatic repulsion and 
can even lead to an attraction. Recently, such short-ranged attractive colloidal systems 
have received renewed attention in terms of their dynamic properties (Nicolai and Cocard, 

2001) . When the strength of short-ranged attraction becomes significantly greater than 
the thermal energy, the system can form a colloidal gel, in which the particle motions 
will be completely jammed even at very low bulk volume fractions. The system behaves 
nonergodic prior to complete blocking of particle dynamics. Therefore, by varying the 
strength and range of the attractive interactions, a reentrant transition of the liquid glass 
line and a glass-glass transition can be realized (Dawson et al., 2001). 

In this paper, we present mainly experimental results obtained with suspensions of 
Laponite at 1.5% wt, / = 7.10"^ M. Very similar results were obtained with suspensions 
of Laponite at / = 10""^ M. For the ionic strength considered here, static and dynamic 
light scattering experiments show that the suspensions form colloidal glasses (Bonn et al., 
1999). 

3 Nonlinear rheology 

Shear experiments were performed on a controlled stress rheometer that can also operate 
in controlled strain mode (Rheologica StressTech), using a Couette geometry or a vane 
geometry with 1 mm gap in the range 50 < 7 < 500 s~^, and a Couette geometry with 1 / 
8 mm gap in the range 500 < 7 < 2000 s~^. The tests were carried out at a temperature 
of 20.0 ± 0.5°C. Evaporation of water or CO2 contamination of the sample was completely 
avoided by covering the sample with vaseline oil and a plate made of Plexiglas. 

At tw = 0, the sample was exposed to a continuous shear strain with a constant rate 
7 between 50 < 7 < 2000 s"^. The viscosity of the suspension increased slowly with 
time and reached a stationary value after some time, as shown in Figure [2l A variation of 
the viscosity of less than 1% per hour was choosen as a criterion for the stationary state 
to be reached, as this is of the order of magnitude of our measurements precision. This 
stationary state was obtained after very long times. To give an order of magnitude, for 
a 4% wt suspension at ionic strengh / = 10""^ M, it was reached in about 50 hours. By 
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increasing the particle concentration from 1.5% to 5%, we observed a slight decrease in 
the time to reach the stationary state. Upon increasing the ionic strength, we were able 
to reach the stationary state of the system faster, as shown in Figure [2l for a suspension at 
1.5% wt, with two different ionic strengths. 

In order to compare the driven aging dynamics to the spontaneous one, we measured 
the complex viscosity modulus of the suspension under shear. The experimental procedure 
was the following : every 30 minutes, we stopped shearing for a very short time (typically 
30 s) and performed oscillatory shear experiments, at a frequency of 1 Hz and at an 
imposed stress of 0.15 Pa, leading to the elastic and loss moduli. This procedure allowed 
us no to perturb the shearing process. The complex viscosity moduli both for the system 
that evolves spontaneously and for the system under shear are shown in Figure [HI As 
the complex viscosity modulus can be taken as a measure of the age of the system, it 
follows immediately that the aging dynamics under flow is very slow compared to the 
"spontaneous" aging dynamics : the complex viscosity modulus increase is slower when the 
system is submitted to an external drive. If the drive was suppressed once the stationary 
state was attained, the "spontaneous" aging process again took place as was observed by 
performing oscillatory shear experiments : the complex viscosity modulus i]* changed very 
rapidly, by 2 orders of magnitude over about 1 hour, in a similar way to what was observed 
for the spontaneous aging process. 

Starting from different aging times t^, a continuous shear strain with the same constant 
rate was applied to the system : the same non-equilibrium stationary state is attained, 
characterized by very similar viscosities. Figure 4 shows the rejuvenation and driven aging 
process under a shear rate 7 = 500 s~\ starting from different aging times. The steady- 
state viscosity thus appears to only depend on the applied shear rate. As the complex 
viscosity can be taken as the measure of the age of the system, the corresponding age of 
the stationary state appears to be power-dependent. When applying the external drive 
to a suspension aged for a long enough time t^, the viscosity decreased in time until the 
system reached its steady state. During this shearing experiment, we measured the complex 
viscosity modulus using the same procedure as explained above. The complex viscosity 
modulus and the dynamic viscosity are shown to evolve in the same direction, they both 
decrease : so-called rejuvenation of the suspension was observed. It is worth noting that, 
even by shearing strongly the suspension during hours, in the range 50 < 7 < 2000 s"^, 
the initial "liquid" state, obtained at = as described in Section II, was never reached 
again, the final state being only determined by the rate of the applied shear strain. 

The stationary viscosities obtained for different values of the applied shear rate from 

= are shown in Figure El The larger the applied shear rate, the smaller the value of 
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Figure 3: Comparison between the aging dynamics under shear (circle) and the spontaneous 
one (square) for a Laponite suspension at 1.5% wt, / = 7.10"^ M: the measured complex 
viscosity modulus varies slowly when the system is submitted to an external drive (7 = 
500 s~^) compared to spontaneous aging. 
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Figure 4: Aging dynamics under a shear strain with a constant rate 7 = 500 s"^ in a 
Laponite suspension (1.5% wt, / = 7.10"^ M). Starting from different aging times — 
(triangle), tu, — 20 min (cross), tw — 40 min (circle), the glass exhibits so-called rejuvena- 
tion (circle) and driven aging dynamics (cross and triangle). The stationary state reached 
only depends on the rate of the applied shear strain. 
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Figure 5: Viscosity as a function of the aging time for a Laponite suspension at 1.5% 
wt, / = 7.10"^ M, for various applied shear rates 50 s~\ 100 s~^ and 500 s~^. The larger 
the applied shear rate, the smaller the value of the stationary viscosity and the faster to 
reach it. 

the stationary viscosity and the faster we reach it. If we consider the typical time necessary 
to obtain the steady state as the time when the viscosity varies less than 1%, we find that 
for the constant respective shear rates 50 s"\ 100 s~^ and 500 s~^, times are respectively 
10 h, 5 h and 1 h. As a result, it roughly scales with I/7. By reporting the values of the 
stationary viscosity as a function of the applied shear rate, a shear-thinning behaviour was 
observed as shown in Figure El It can be accurately characterized by a power law 77 oc 7"", 
with a = 0.60 ± 0.10 for the different suspensions of Laponite considered here. 
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Figure 6: Stationary viscosity as a function of the applied shear rate for a suspension 
at 1.5% wt, / = 7.10"^ M (square) and a suspension at 3.7% wt, / = 10~^ M (circle). 
We observe a shear-thinning behaviour characterized by a a power law rj oc 7~", with 
a = 0.60 ± 0.10. 
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4 Discussion 



From a theoretical point of view, a number of models have been introduced recently to 
account for the interplay between the drive and the relaxation of the system. Sollich 
and co-wokers (1997) extended Bouchaud 's trap model (Bouchaud, 1992; Monthus and 
Bouchaud, 1996) to driven systems, leading to the 'Soft Glassy Rheology' (SGR) model 
(Sollich, 1998; Fielding et al, 2000). Other models that capture many of the observed rheo- 
logical properties have been introduced (Hebraud et al, 1998; Derec et al., 2001, Coussot et 
al., 2002). Cugiiandolo et al. (1997) suggested an approach to non-equilibrium systems, in 
which the non-equilibrium state was generated by 'stirring' the system. Further theoretical 
studies on driven mean-field disordered systems and numerical studies on Lennard- Jones 
glasses allowed for a detailed description of the non-equilibrium behaviour, including rheo- 
logical properties (Berthier et al, 2000; Barrat and Berthier, 2000). All these models were 
shown to account well for a number of rheological behaviors in soft glassy materials. In 
particular, as a result of these theoretical predictions, it is expected that if we inject power 
into aging systems, there is the possibility of stabilizing the age of the system in a power 
dependent level : the younger the larger the power input (Bouchaud et al., 1996; Horner, 
1996; Sollich, 1998; Kurchan, 1999; Berthier et al, 2000; Derec et al., 2001). The power 
injection into the system derives from the applied shear rate. For a vanishingly small drive, 
however, the slow mode relaxation time, assumed to be proportional to the viscosity, again 
diverges (Berthier et al., 2000). 

These theoretical predictions are confirmed by the experimental results described above. 
The driving force acting on the system results in the suppression of the aging process. The 
larger the shear rate, the smaller the stationary viscosity that means the younger the glassy 
system, as the complex viscosity is known to increase with time t^, and can be therefore 
taken as a measure of its age. Once the drive is suppressed, the spontaneous aging process 
again takes place. From our experimental results, it appears that the external drive leads 
to a drastic slowing down of the spontaneous aging dynamics or, in some cases, in the 
rejuvenation of the system. The occurence of rejuvenation or "driven" aging depends on 
both the external drive power and the age of the system when applying the external 
drive. This is completely consistent with phenomenological models in which the final 
stationary state results from the competition between the spontaneous aging dynamics - 
that strenghtens the interactions between particles and leads to less and less accessible 
configurations in phase space as time evolves - and the mechanically induced rejuvenation 
of the configurations - process going against the spontaneous aging and leading to its 
slowing down. Coussot and co-workers (2002) have shown that, under imposed stress a 
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bifurcation in rheological behavior occurs: for small stresses, the viscosity increases in time: 
the material eventually stops flowing. For slightly larger stresses, the viscosity decreases 
continuously in time: the flow accelerates. Thus, the viscosity jumps discontinuously to 
infinity at the critical stress, implying also there is a range of shear rates for which the 
fiow is unstable. Here, we work at sufficiently high imposed shear rates so that the fiows 
are stable; no shear banding is observed. 

In the experiments, the typical time to reach the stationary state is surprisingly long. 
Our experiments show that this typical time scales with I/7, result which confirms theoret- 
ical expectations and numerical simulations (Berthier et al, 2001). However, the prefactor 
of 1/7 is of the order of 10^ in a typical experiment and suggests that the applied shear 
disrupts very large scale structures. We observed that this typical time is of the same order 
of magnitude than the time necessary for the spontaneous aging dynamics to slow down 
dramatically. This fact strengthens the idea of a competition, between spontaneous aging 
and the external drive. Furthermore, these steady-shear experiments disprove a common 
idea on Laponite suspensions : when taking an "old" suspension, it is not sufficient to stir 
the system during a few minutes in order to obtain the initial state at tu, = 0, for the shear 
strain rates considered here. 

The shear rejuvenation was studied by using Diffusive Wave Spectroscopy (DWS) in a 
suspension of Laponite (Bonn et al., 2002). The shear rejuvenation was shown to have a 
very large effect on the microscopic dynamics. Upon increasing the shear rate, the steady- 
state relaxation time decreases. In the experiments reported in this paper, the complex 
viscosity modulus associated to the stationary state was also shown to decrease with the 
applied shear rate. This is consistent with the DWS experiments, where the measurements 
allowed us to relate the microscopic diffusion dynamics to the macroscopic viscosity of the 
system. 

In the SGR model, the rheological response of the system depends on the distance from 
the glass transition temperature x. If we assume that we are in the immediate vicinity 
of the glass transition, x ~ 1, the shear-thinning exponent found experimentally is in 
reasonable agreement with the theory. Berthier et al. (2000) studied, within the framework 
of the mean-field (or mode-coupling) approximation, the infiuence of an external drive on 
a system undergoing a glass transition. Above the glass transition, for high enough shear 
rates, the slow mode relaxation time ta of the system scales as ta oc 7"^/^. As the dynamic 
viscosity 77 is assumed to be proportional to the slow mode relaxation time, the system has 
therefore a shear-thinning behavior characterized by the power law rj oc 7"^/^. Very similar 
results are obtained from molecular dynamics simulations, where the viscosity is directly 
extracted from the applied shear rate (Barrat and Berthier, 2000; Yamamoto and Onuki, 
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1998), with T) oc 7"°-^ and 77 oc 7"°-^ respectively. The shear-thinning behavior observed 
experimentally (an exponent between —0.5 and —0.7) is in fair agreement with the above 
theoretical and numerical predictions. 

In conclusion, we have studied experimentally the nonlinear rheological behavior of 
a glassy system under an external drive. Qualitatively, our experimental results confirm 
all the findings predicted by the glassy models. They show the competition between the 
external drive, leading to a mechanically induced rejuvenation of the configurations, and 
the aging dynamics, leading to a slowing down. Depending on the competitive antagonists, 
driven aging or rejuvenation may be observed. In addition, the shear-thinning viscosity 
observed in our experiments agrees quantitatively with at least one simulation result. Since 
the value of the exponent depends on the exact model that is chosen, the differences with 
the model calculations may not be very important. The key result is that non-equilibrium 
models for glassy systems are able to predict a non-Newtonian behavior without taking 
into account the specific interactions between particles. 
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